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1 Introduction 

Localisations of intracranial electromagnetic sources 
with magnetoencephalography (MEG) are usually 
computed with mathematical models that 
approximate the skull volume by means of a 
homogenous sphere fitted to the surface of the skull 
where the sensors record strong magnetic field 
gradients. 

In a study of the validity of the spherical model, 
localisation accuracy was determined using an 
artificial dipole in various regions of head phantoms 
in different dipole directions [1]. Localisation errors 
were found to be small and independent of dipole 
direction when dipoles were placed under spherical 
surfaces, e.g. in parietal regions. In contrast, when 
dipoles were placed under non-spherical surfaces, 
e.g. in frontal and temporal regions, errors were 
large or localisation not possible due to poor 
correlation between observed and theoretical fields 
for many dipole directions. 

However, when the recordings from the sensors 
contralateral were added and localisations computed 
with the skull volumes approximated by means of a 
boundary element model (BEM), these errors were 
largely diminished [2]. 

To determine whether the improved localisation 
accuracy was due to the more realistic description of 
the skull volumes with the BEM or to the use of 
bilateral recordings, localisations based on bilateral 
recordings using the spherical model and unilateral 
recordings using the BEM were compared with 
previous localisation results. 

2 Material 

A plastic phantom with a closed internal surface 
similar to the internal surface of a real skull was 
moulded using a form made of wooden plates 
modelled from 3 mm successive coronal MRI slices 
of a human head. It was filled with a saline solution 
with homogenous conductivity. The phantom was 
equipped with a ball-joint through which a plastic 
tube carrying a dipole was inserted, to allow the 


dipole to be positioned in different regions of the 
volume (see Fig. 1). 



Figure 1: Artificial dipole set-up for recordings of 
generated magnetic fields. The dipole was inserted 
in a skull shaped phantom filled with a 0.9% saline 
solution and placed between the upper and lower 
dewar of a 74-channel biomagnetometer system. The 
positions of the dipole instrument and the phantom 
were kept fixed by means of a Delrin construction. 
The dipole orientations were set with the dials at the 
outer end of the tube without altering the centre of 
the dipole. The dipole could reach all regions of the 
phantom by directing and sliding the tube through 
the ball-and-socket-joint; in the example shown, it is 
placed in the left temporal region. 

The dipole was constructed from two carbon rods 
symmetrically inserted into a plastic cylinder 
leaving the surfaces of the outer tips exposed. Each 
of the exposed areas was 3 mm 2 and the distance 
between them was 14 mm. A pair of twisted copper 
leads was used to connect the carbon rods to a 
current source through the plastic tube. 


The dipole was mounted on the tip of the plastic 
tube in a fixture allowing graded changes of the 
dipole’s direction without altering the position of the 
dipole’s centre by means of two dials at the other 
end of the tube outside the phantom. The current 
source, ~ 0.1 mA, was controlled by a 10 Hz 
sinusoidal voltage and the magnetic fields measured 
with a dual 37 channel MEG system (Magnes II, 
BTi) where one of the sensors was positioned to 
record bipolar magnetic fields with strong gradients 
while the other was placed on the contralateral side. 

3 Methods 

The magnetic fields generated by the dipole were 
measured with the dipole in different regions of the 
skull shaped phantom. The regions were comparable 
to frontal, mesial temporal and parietal regions in 
the corresponding real head. 

For each position in the phantom, the dipole was set 
in 103 different directions evenly distributed within 
a sphere and the magnetic fields were measured at 
each dipole direction. 

The measured signals were filtered using a bandpass 
filter with cut-off frequencies of 5 and 15 Hz. The 
DC component of the filtered signals was then 
subtracted removing any residual offset. 

For each measured field the dipole localisation was 
calculated using 1) software supplied by the 
manufacturer of the MEG system based upon the 
spherical model, capable of utilising unilateral 
recordings only [3] and 2) software based upon a 
boundary element model adapted to unilateral as 
well as bilateral recordings [4]. The latter software 
was modified to also include a spherical model 
capable of utilising bilateral recordings. 

For all four methods, starting points for the iterative 
search algorithm were determined by calculating 
localisations at grid points with 10 mm spacing 
covering the whole region of interest and selecting 
the grid point with the best fit. 

Focalisation results were accepted if the goodness- 
of-fit between the field from a theoretical source and 
the measured field were >0.98. 

4 Results 

In Figs. 2 and 3 the results are illustrated by 
comparisons between localisations based on 
spherical and boundary element models with uni- 
and bilateral recordings in frontal and mesial 
temporal regions. 

In the frontal region when localising using the 
spherical model with unilateral recordings, the 
localisations were scattered between the centre of 


the approximating sphere and the volume surface 
close to the recording sensors (see Fig. 2). 


Frontal region 



Figure 2: Localisation accuracy in frontal region 
using spherical and boundary element models with 
uni- and bilateral recordings. The phantom is shown 
in sagittal and coronal projections and true dipole 
position is indicated by the intersection between co¬ 
ordinate axes with tickmarks at 1 cm intervals. 
Magnetic fields were generated with the dipole's 
centre in fixed position while it was turned into 103 
evenly distributed directions. Trials to localise the 
dipole from these fields were made. Accepted 
localisations are given by red dots. For spherical 
models, green crosses indicate centre of the 
approximating sphere. 



















When localising with the BEM and unilateral 
recordings some localisations also gave relatively 
large errors. Localising with bilateral recordings, 
however, gave small errors regardless of the model, 
although errors were less with the BEM. 

In the mesial temporal region when localising using 
the spherical model with unilateral recordings, the 
localisations were also scattered but without any 
apparent pattern (see Fig. 3). 


Mesial temporal region 
Sagittal Coronal 




Figure 3: Localisation accuracy in mesial temporal 
region using spherical and boundary element 
models with uni- and bilateral recordings. Figure 
details as in Fig. 2. 

Similar results were obtained when localising with 
the BEM using unilateral recordings and with the 


spherical model using bilateral recordings. Thus, in 
this region only when localising with the BEM and 
bilateral recordings were the errors small. 

In Fig. 4 graphs show the distribution patterns of 
localisation errors >10 mm when using the four 
different localising methods in frontal, mesial 
temporal and parietal regions, respectively. It also 
quantifies the number of localisations with errors 
>10 mm and <10 mm as well as the number of 
localisations rejected due to goodness-of-fit < 0.98. 
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Figure 4: Distribution of localisation errors > 10 
mm using spherical model (Sph.) and boundary 
element model (BEM) with unilateral (U) and 
bilateral (B) recordings. Below the graphs, number 
of accepted localisations with errors >10 mm and 
<10 mm as well as number of localisations rejected 
due to goodness-of-fit < 0.98 are given. 

With the bilateral BEM most localisations (94 - 96 
%) were accepted and only few (0-2 %) had errors 
larger than 10 mm (max. 14 mm). With the bilateral 
spherical model many localisations were rejected. In 
remaining localisations, the localisation errors were 
often large (> 10 mm) but rarely very large (>20 
mm) as with the unilateral spherical model. The 
unilateral BEM gave the opposite result, with most 
localisations accepted. Relatively few of these (6 - 
23 %) had errors larger than 10 mm but among 
these, very large errors (20 - 40 mm) were common. 



































5 Discussion 

In frontal and temporal regions where the head 
surface is non-spherical, the spherical model with 
both uni- and bilateral recordings gave localisations 
that were not acceptable or had large errors in many 
of the dipole directions tested. In contrast, when 
using the BEM with bilateral recordings, the dipole 
could be localised with small errors in nearly all 
dipole directions. Satisfactory accuracy, however, 
was not obtained when using the BEM with 
unilateral recordings. In the parietal region, where 
the head surface is close to spherical, errors were 
small except for the BEM with unilateral recordings. 
Previous reports on localisation accuracy of the 
spherical model have usually reported smaller errors 
than those reported in the present study. Some of the 
reports have also provided real measurements but 
only with relatively few dipole positions and 
directions [e.g. 5-7]. Because brain malformations 
and lesions in the brain of epileptic patients may 
reduce the possibility of predicting source 
directions, the accuracy of all dipole directions 
should be taken into consideration. Other reports 
have used mathematical simulations instead of real 
measurements as a mean of determining the 
localisation accuracy [e.g. 8-10]. While it is less 
difficult to test the various positions and directions 
with simulations, some of the complexity is likely 
lost, thus effecting the accuracy of determinations. 

In two recent reports, Leahy et al. [7] and Crouzeix 
et al. [10], localisation accuracy with the BEM and 
the spherical model were compared. Although the 
absolute differences between models were less than 
found in the present study, the BEM was determined 
to be superior in both studies. Since Leahy et al. 
concentrated on sources in visual and somatosensory 
areas, i.e. in regions close to spherical skull surfaces, 
it is reasonable to also find a relatively high 
accuracy with spherical models. As expected, the 
two dipoles in other regions had larger errors. 
Although Crouzeix et al. considered dipoles in most 
regions with various directions they elected to use 
the exact positions of the original dipoles as starting 
points for their iterative search algorithm. This may 
create bias to end searches at a local minima close to 
the starting points and might well explain the rather 
small difference observed between using the BEM 
and the spherical model. 

The results of this study strongly suggest that 
bilateral recordings of the magnetic fields and the 
use of a realistic volume model such as the BEM are 
needed to ensure a high source localisation accuracy 
with MEG in all regions of the brain. 
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